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The dynamical response to a pulsed magnetic field has 
been studied here both using Monte Carlo simulation and 
by solving numerically the meanfield dynamical equation of 
motion for the Ising model. The ratio R v of the response 
magnetisation half-width to the width of the external field 
pulse has been observed to diverge and pulse susceptibility 
Xp (ratio of the response magnetisation peak height and the 
pulse height) gives a peak near the order-disorder transition 
temperature T c (for the unperturbed system). The Monte 
Carlo results for Ising system on square lattice show that R p 
diverges at T c , with the exponent vz = 2.0, while \v shows 
a peak at T^, which is a function of the field pulse width 
St. A finite size (in time) scaling analysis shows that = 
T c + C(8t)~ 1/x , with x = vz ^ 2.0. The meanfield results 
show that both the divergence of R and the peak in \v occur 
at the meanfield transition temperature, while the peak height 
in Xp ~ (8t) y ', y — 1 for small values of St. These results also 
compare well with an approximate analytical solution of the 
meanfield equation of motion. 

PACS number(s): 05.50.+q 



I. INTRODUCTION 

The dynamic response of the Ising systems has recently 
been studied extensively employing computer simulations 
[1]. In particular, the study of dynamical response of 
Ising systems to oscillating magnetic field [2,3] has led 
to many intriguing dynamic phenomena, like dynamic 
hysteresis and the fluctuation induced dynamic symme- 
try breaking transitions in (low e.g., one, two or three 
dimensional) Ising system in presence of an oscillating 
field. Acharyya and Chakrabarti also noted [3,4] some 
anomalous behaviour in the growth of pulse susceptibil- 
ity, in Ising systems under pulsed magnetic fields of finite 
durations. 

Usually, when a cooperatively interacting thermody- 
namic system in equilibrium is perturbed (with the per- 
turbation having a step function like variation with time), 
then the relaxation of the system (to the equilibrium ap- 
propiate to the perturbed state) is observed to follow 
the common Debye type form with a single relaxation 
time. The standard (Debye) form for any response func- 
tion (say magnetisation of Ising system) m(t) is 



equilibrium value, A is a constant. As the critical tem- 
perature is approached r shows a critical slowing down; 
r diverges at the critical temperature T c : 



~ (T-T c y 



(1.2) 



where £ is the correlation length, z is the dynamic expo- 
nent and v is the correlation length exponent [1]. 

Here, we have investigated in details the response of 
pure Ising systems to pulsed magnetic fields of finite du- 
ration, using Monte Carlo simulations for two dimen- 
sional Ising systems, and solving numerically the mean- 
field equation of motion. We have studied the response 
behaviour for 'positive' pulses, where the pulsed field is 
in the direction of the sponteneous magnetisation in the 
ordered phase. In the disordered phase, of course, this 
notion is immaterial. One can also study the effect of 
'negative' pulses on the sponteneous order, where the 
field direction is opposite to sponteneous magnetisation. 
Although many intriguing features of the domain growth 
etc are expected for such 'negative' pulse problem, we 
confine here to the study of 'positive' pulses only. We 
have measured the ratio R p of the response magnetisa- 
tion (pulse) half-width (At) to that (St) of the pulsed 
external field, and the ratio Xp °f the response magneti- 
sation peak height (m p ) to the field pulse height (h p ) giv- 
ing the pulse susceptibility. The temperature variation 
of these two quantities for various pulse width durations 
and heights of the external field, have been investigated. 

We find that for weak pulses, while the width-ratio 
R p diverges at the order-disorder transition point T c of 
the unperturbed Ising system, the pulse susceptibility Xp 
does not diverge and for low dimensions, e.g., in two di- 
mension studied here, a smeared peak in \p occurs at an 
effective T c e , which approaches T c as the field pulse width 
St increases (x P diverges at T c as St —> oo). A meanfield 
analysis for the absence of divergence of Xp f° r finite St 
has been developed. Also, a finite time scaling analysis, 
similar to Fisher's finite size (length) scaling [5], has been 
developed and compared with the observation of the ef- 
fective transition temperature T c e with the external field 
pulse width St. 

We have organised this paper as follows: In section 
II, the model and the simulation techniques have been 
described. In section III, the results are given. The paper 
ends with concluding remarks in section IV. 



m(t) ~ m(oo) + A exp(— t/r), 



(1.1) 



where r is the relaxation time, m(oo) denotes the new 
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II. THE MODEL AND SIMULATION 

A ferromagnetically interacting (nearest neighbour) 
Ising system in presence of a time dependent magnetic 
field can be described by the Hamiltonian 



meanfield dynamical equation of motion for the kinetic 
Ising system: 



(2.1) 



where sf's are spin variable having their value ±1 and 
h{t) is the time varying longitudinal magnetic field. Here, 
we have considered the time variation of h(t) as follows: 



h(t) = h p for t < t < t 
— elsewhere, 



St 



(2.2) 



where h p is the amplitude of the field and St is the dura- 
tion or the active period of the external field. 

In our simulation, we have considered a 500x500 
square lattice in two dimension. At each site of the lattice 
there is a spin variable sf = (±1). We update the lat- 
tice by stepping sequentially over it following the Glauber 
single spin flip dynamics. One such full scan over the lat- 
tice is unit time step (Monte Carlo step or MCS) here. 
First, we allowed the system to reach the equilibrium (at 
any temperature T) and only after that the magnetic field 
h(t) has been switched on (t is thus much larger than the 
relaxation time of the system) . One can also use random 
updating sequences. However, it takes more time (MCS) 
to stabilise the system. We expect all the features of the 
response studied here remains qualitatively unchanged 
for random updating sequences, and we give the results 
here for sequential updating only. We have measured the 
maximum height (above the equilibrium value) m p and 
half-width At of the response magnetisation. Here, as 
mentioned before, the following two quantities have been 
defined to characterise the response of the system: 



pulse width ratio R p 



At 
~St 



and 



•1 .T " l V 

pulse susceptibility Xp = -rf-- 

tin 



It may be noted that Xp reduces to normal static suscep- 
tibility as hp — > and St — > oo. At each temperature for 
fixed hp and St, the numerical values of R p and Xp are 
obtained averaging over 20 random Monte Carlo reali- 
sations (initial seed). We have studied the temperature 
variation of these two quantities. The results are given 
and discussed in the next section. The observation for 
finite peak in Xp a ^ an effective transition temperature 
Xj? (which converges to T c as St — > oo) is analysed 
in view of finite time scaling behaviour, as mentioned 
before. 

To study the similar response in the case where the 
fluctuations are absent, we have considered the following 



dm ( m + hit) 

rn = — m + tanh 

dt \ T 



(2.3) 



Here To is the microscopic relaxation time, m is the aver- 
age magnetisation (in the meanfield approximation), h(t) 
is the time varying pulse fi eld having the same time varia- 
tion described in equation 2A and T denotes the temper- 
ature. We have solved numerically the above equations 
using fourth order Runge-Kutta method. We have eval- 
uated the pulse width ratio R p and pulse susceptibility 
X P at various temperatures for fixed pulse width St and 
height hp. The numerical results are given in the next 
section, where we have compared the results for finite 
peak height in Xp (at T c = 1), with an approximate an- 
alytic estimate of Xp m such cases. 



III. RESULTS 
A. Monte Carlo studies 

As mentioned before, our results here are for two di- 
mensional Ising system on square lattice of size 500x500. 
We applied an external field of amplitude h p for a du- 
ration St after bringing the system into a steady state. 
For this, the typical number of Monte Carlo steps re- 
quired for this size of the lattice chosen here, is observed 
to be of the order of 10 6 . The response magnetisation 
has amplitude m p (measured from the equilibrium value) 
and a half- width At. Fig. 1 shows a typical time vari- 
ation of magnetic field h(t) and the corresponding re- 
sponse magnetisation m[t). The dynamical response is 
characterised by two quantities; the width-ratio R p and 
the pulse-susceptibility Xp- The temperature variation of 
these two quantities has been studied. Fig. 2 shows the 
temperature variations of R p and Xp f° r fixed values of 
h p (= 0.5) and for three values of St (= 5, 10, 25 MCS). 

Since m p is bounded from above, a large h p would 
saturate m p and hence Xp becomes small (due to the 
saturation). Also, for extremely small h p , it becomes 
difficult to identify m p from the noise, and hence the 
estimate of Xp becomes erroneous. We found h p = 0.5 
to be well within the above optimal range. From the 
figure it is clear that R p has a sharp divergence at Tr(= 
2.30, somewhat larger than T c , the Onsager transition 
temperature, due to the small size of the system) almost 
irrespective of the values of the pulse width St. But Xp 
shows peak at different points (significantly above T c ~ 
2.27) depending upon the values of the field pulse width 
St. As St value increases, it is observed that the peak 
shifts towards T c from above (and also the peak height 
grows) . 

Let us try to understand why the width-ratio R p di- 
verges at T c , while height-ratio or pulse susceptibility Xp 
shows peak at some higher value T^{8t) depending upon 
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the value of St: T£(St) — > T c as St — » oo. The pulse-like 
perturbation probes the response of the system at finite 
frequencies. Consequently, the x P (= m p /h p ) can not di- 
verge as the response magnetisation height m p is not an 
equilibrium value corresponding to the pulse height h p ; 
rather the m p results are bounded by the time window 
of width St. On the other hand, the response will take 
its own relaxation time to come to its equilibrium value 
(irrespective of the value of St), when the field is switched 
off (at to + St). This leads to the divergence of At, due 
to critical slowing down, as T approaches T c . 

The sharp divergence of width-ratio R p is identified as 
the consequence of critical slowing down and the point 
of divergence is the critical temperature T c for the ferro- 
para transition. In fact, since the relaxation after the 
withdrawal of the pulse is unrestricted by the pulse- 
width, we can assume that At ~ r ~ |T — T c \~ vz . We 

therefore plot in Fig. 3, R p 1 versus T and find a 
straight line plot with vz = 2.0 in two dimensional case. 
This compares well with the previous estimates of the 
value of vz [6] . 

As the growth of the height of magnetisation response 
(and its maximum value m p ) is very much bounded by 
the time window St of the applied field pulse, the anoma- 
lous behaviour of Xp (having finite peak at a shifted tem- 
perature Tf!(St)) may be considered to be due to the finite 
size (in time) effect. Similar to the finite size (in length) 
scaling theory of Fisher [5], where a finite size system 
shows effective (non-singular or non-divergent) pseudo- 
critical behaviour at Tf;(L) when the correlation length 
£ becomes of the order of the system size L, we suggest 
a finite time (St) scaling behaviour here for Xp'- K the 
relaxation time r ~ £ z ~ \T — T c \~ vz , where v is the 
correlation length exponent and z is the dynamic expo- 
nent, then Xp would show peak at the temperature T c e 
here when r(T c e ) - St or \T*(5t) - T c \~ vz ~ St, or 



T*(5t) ~ T c + C(St)- 1 ' vz 



(3.1) 



where C is some constant. In fact, Fig. 4 shows that the 
effective peak position T c e indeed approaches T c as St — » 
oo. The inset shows the plot of T% with (St) -1 ' 3 *, which 
gives a straight line for x = vz = 2.0. This again suggests 
vz = 2.0 and also the extrapolated value of T c becomes 
about 2.29, which compares well with the Onsager value, 
comparable to the previous estimate. 



B. Meanfield results 

We have sol ved the meanfield equation for response 

fourth order Runge-Kutta 



magnetisation 
method. Fig. 



2.3 



using 



1 shows the typical variation of response 
magnetisation and field. Here also we have measured the 
width-ratio and the pulse susceptibility and studied the 
temperature variation of these two quantities. 

Fig. 5 shows the temperature variation of R p and Xp 
for different values of St. R p diverges and Xp peaks at 



the same order-disorder transition point (T c = 1 here). 
We have also studied the variation of the maximum value 
of Xp (x™^ at T = T c ) with respect to the duration of 
the pulsed field. It may also be mentioned that the peak 
height was found to increase with increasing pulse width 
St, and Xp diverges as St — ► oo: x r p ax ~ (8t) v ;y = 1.0 
Fig. 6 shows the variation of x™ a2: with St. 

In order to comprehend these observations, we solve 
the equation 2.3 in a linearised limit (large T and small 
h p ; specifically T > 1, h p — > 0). In such limit, the equa- 
tion of motion becomes 



r ^ = -em + h(t)/T; e=(T-l)/T. (3.2) 

One can solve the above equation using a form m(t) = 
moe _ */ T which gives 



at 



giving t/tq = e , and 



T 

-1 



T ' 



dmo e -t/r = 
dt 



= h(t)/T. 



(3.3) 



(3.4) 



Integrating the last equation for h(t) = h p for a finite 
time width St and h(t) = elsewhere, one gets m p ~ 
h p St/(T c T ) at T = T c = 1 (when r -> oo). This gives 

X p aX = X P (Tc) = m p (T c )/h p ~ St/(T C T ) ~ St/r Q . (3.5) 

We have checked the above linear relationship of X™ ax 
with St for extremely small values of pulse field amplitude 
h p (see Fig. 6). 



IV. SUMMARY 

We have studied the Glauber (order parameter noncon- 
serving) dynamics of an Ising system under a time vary- 
ing external magnetic field, when the field is applied as 
a pulse of finite time width, after the system reaches the 
equilibrium. The time variation of the response magneti- 
sation is studied as a function of pulse width St, height 
hp and temperature T of the system. We have measured 
specifically R p = At /St and Xp = Tn p /h p , where At is the 
time width of the response magnetisation and m p is the 
maximum height of the response magnetisation above its 
equilibrium value. 

Our computer simulation results for square lattice 



showed R n 



T 



with T c = 2.30, the On- 



sagar value, and Xp nas a peak x!p ax at T£ > T c , 
such that a finite size (in time) scaling behaviour is ob- 
served: T c e = T C + C(8t)- 1 / vz , with vz 2.0. The nu- 



merical solutions of the mean field equation 2.3 showed 
Rp ~ 1/(T — 1) and the pulse susceptibility peak value 
^max ^ occuring at T = T c = 1. Theoretical analysis 
for the finite size (in time) scaling behaviour (of in 
the Monte Carlo case) and for the peak x™ aa; (<5t) (in the 
mean field case) are also given. 
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Figure Captions 

Fig.l. Time variation of magnetic field (h(t)) and the 
response magnetisation (m(t)) in the Monte Carlo case. 
For hp = 0.5 and St = 50. 

Fig. 2. Temperature variations of (a) R p and (b) \p f° r 
different values of St in the Monte Carlo case. The symbol 
circle is for St = 5; the square is for St = 10 and the cross 
is for St = 25. 

Fig. 3. Variation of Rp 1 ^ vz against T, with vz = 2. 
Fig. 4 Variation of T c e with respect to 1/St in the Monte 
Carlo case. Inset shows the variation of Tj? with respect 
to (5t)-^ x ; x = 2. 

Fig. 5 Temperature variation of (a) R p and (b) \v f° r 
different values of St in the mean field case. The symbol 
triangle is for St = 8; the plus is for St = 16 and the cross 
is for St = 32. 

Fig. 6. Variation of Xp lax with respect to St in the mean 
field case. 
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